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Deep-inelastic lepton-hadron scattering is analyzed in asymptotically free gauge theories 
of the strong interactions. The renormalization-group equations for the coefficients of the 
twist-two operators in the Wilson expansion are reviewed. A careful treatment of the mix- 
ing of operators with identical quantum numbers and dimensions is given. The relevant 
anomalous dimensions of the twist-two operators are calculated to second order in perturba- 
tion theory. These are used to calculate the asymptotic q? behavior of the moments of the 
structure functions. It is shown that the approach to the asymptotic limit is logarithmic, 
that Bjorken scaling is violated by powers of In(—q’), and that the naive light-cone or parton- 
model relations for the moments of the structure functions are true asymptotic theorems. 

A new sum rule for the first moment of F,, in terms of the energy-momentum tensor, is 
derived. An example of a function whose moments have roughly the correct asymptotic 
q? behavior is constructed. The q? behavior of the structure functions for a given x is 


discussed. 


I. INTRODUCTION 


In a recent paper’ we have constructed a class of 
gauge theories of the strong interactions, which 
have the remarkable feature of being “asymptot- 
ically free.” The primary motivation for this 
proposal is the evidence that Bjorken scaling re- 
quires an asymptotically free theory,’ that only 
non-Abelian gauge theories can be asymptotically 
free,* and that indeed many non-Abelian gauge 
theories are asymptotically free.*~* Deep-inelastic 
scattering is therefore the natural arena in which 
to test our theories. In this paper we shall dis- 
cuss in detail the properties of lepton-hadron 
scattering in asymptotically free gauge theories of 
the strong interactions. 

This paper is a sequel to Ref. 1 and should be 
read in conjunction with it, although the phenom- 
enological discussion of Sec. III can be understood 
independently (hereafter Ref. 1 will be referred 
to as paper I, and the prefix I refers to equations 
of Ref. 1). The general features of deep-inelastic 
scattering in asymptotically free theories were 
already described in paper I. These include the 
logarithmic approach to scaling, the calculable 
logarithmic deviations from Bjorken scaling, and 
the validity of the naive or light-cone parton-model 
relations. 

In Sec. II of this paper, we discuss in some de- 
tail the application of renormalization-group tech- 
niques to the Wilson expansion. This analysis has 
appeared in many other places’~* and is included 
here for the sake of completeness. In particular 
we discuss the mixing of operators with the same 
quantum numbers and dimensions. In gauge the- 
ories this mixing is particularly annoying since 


K- 


it would appear that “ghost” operators (i.e., oper- 
ators involving Feynman-Faddeev-Popov ghost 
fields!” ™!) mix together with ordinary operators. 
We argue that this mixing can be ignored. This 
claim is further substantiated by a calculation, 
which appears in Appendix A, performed in a gauge 
which is free of Faddeev-Popov ghosts. 

In Sec. III we calculate the anomalous dimensions 
of the relevant operators in the Wilson expansion.!” 
Some of these were already presented in paper I. 
We derive the asymptotic form for the moments of 
the structure functions, as well as the various 
relations and sum rules satisfied by these moments. 
A sum rule for the first moment of the structure 
functions, the “energy-momentum-—tensor sum 
rule,” is derived. An example of an explicit func- 
tional form for the structure functions, with rough- 
ly the correct asymptotic behavior, is presented 
and the general features of this function are dis- 
cussed, ee: 

Section IV contains some concluding remarks. 


II. THE RENORMALIZATION -GROUP APPROACH 
TO THE WILSON EXPANSION 


In the deep-inelastic scattering of a lepton 
off a hadron one measures the Fourier transform 
of the commutator of electromagnetic or weak 
currents. We define the standard structure func- 
tions as follows: 


ar J 49 eI (a5 9), (05 ~ 4 99] [DD 


-Pube FOD (y, q3) -E rE (y, q3) 


apen (v, q’) tereg (1) 


980 
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where p is the momentum of the hadron, q is the 
momentum transfer to the hadrons, v=p-°gq, spin 
averages have been taken, and we have suppressed 
terms proportional to q, and q, which give terms 
proportional to the lepton mass when contracted 
with the leptonic currents. The label a denotes 
the SU(3) x SU(3) character of the current. The 
Bjorken limit corresponds to 


y/ y? te, 


gees fixed, Q?= o0 2 
“9y ’ pera «ee (2) 


and probes the commutator for lightlike values of y. 

The structure of the product of local operators 
at short, and at lightlike, distances is given by 
Wilson’s operator-product expansion. In the case 
of interest this expansion takes the form 


a\? 1 
J,(a38 9) Jy (05 -49)=48,0( 2] aa 2 > cf (a, b; y?” — icyo) On, oe p,O)ytte ++ yn 


vy a n=0 


1 co 


. ð ; 
+2 ieuvon gyi jaore, Da dos Chs (4s b; Y? = FEV) Op,» 


n=0 


In this expansion we have neglected, as in Eq. 
(1), gradient terms (involving 3/ay" or a/ay’). 
The operator O$. . 1, (x) has spin k (traceless 
and aupametric) vad dimension n +2 (or twist 
= dimension — spin =2). The label z denotes the 


various operators of equal twist which might occur 


in the expansion. C{)}(k=1,2,3) are c-number 
functions of x? and the coupling constants of the 
theory. They have been normalized to be dimen- 
sionless, so that naively (as in free-field theory) 
they would be nonsingular as x?~0. The three 
dots refer to other operators of higher twist, 
whose contribution to the structure functions is 
suppressed, to any finite order of perturbation 
theory, by powers of g?. (In an asymptotically 
free theory this suppression is guaranteed to all 
orders in perturbation theory by the vanishing of 
the anomalous dimensions of all operators at the 
fixed point g=0. One has logarithmic, but not 
power, corrections to naive scaling.) 

The structure functions F, in the deep-inelastic 
region are determined by the behavior of C$} for 
small x, and by the hadronic matrix elements 
of the operators oO! 


(p lO}, "Uy (0) | P? spin average 
asi 
=i"— py * tebu, Mares. (4) 


One cannot simply take the Fourier transform of 
Eq. (3), since after Fourier transforming it does 
not converge in the region of interest (0 <x <1). 
Instead the Fourier transforms of Cr are deter- 
mined by the moments of the structure func- 
tions.°° 

In fact 


> 2, Cfa (a, b; y? = i€) Oy, - 


. p, (0)y"1e ee y"n 


TR a Se E (3) 


1 
f dex" FE” (x, q?)= SD CHP (a, b; qM", 
(0) 
T dx x” F$” (x, q°)= 2 Ci (a, b; q?°)\M?3*?, (5) 
f dx x" F$” (x, q?)= 2 Cis (a, b; q? M3, 
where 
~ : ə n 
H (a, b; q) = 8il (~ 55) 
q 
x fas ye 


In a free-quark model the operators Of... p, 
are simply 


et y HAC biy 2) , (6) 
y? —1EYo 


Oh er u (*)=S P(x) Yu, One’ . ðu A (1 +y) P(x) , 
| (7) 


where S denotes symmetrization of the indices 
Hı°** Hp. In that case the functions C are con- 
stants, independent of q°, which can be deter- 
mined by the light-cone commutator of the cur- 
rents (see Ref. 14). 

In an interacting theory ile functions C™ will 
be nontrivial functions of g? and the coupling con- 
stants of the theory (g). To determine the g? de- 
pendence of Č™ we employ the renormalization- 
group equations. Let us outline the derivation of 
these equations for the Wilson coefficients. Con- 
sider the contribution of an operator O™” (x) to the 
short-distance expansion of the product of A(x) 
and B(x) (we suppress all tensor and quantum- 
number labels): 
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A(x) B(- x) X C” (x?, ORR ep (O) xtian, 


(8) 


If O™ (x) is the dominant operator for short dis- 
tances (i.e., the operator of smallest dimension, 
which we assume for the moment to be unique) 
then the short-distance behavior of C™ can be 
determined by calculations performed as if all 
masses and dimensional coupling constants were 
zero (for details see paper I). In that case the 
only dimensional parameter in the theory is the 
subtraction point u introduced to perform the re- 
quired renormalization. A change in u can be 
reabsorbed by a change in the coupling constant 
and the scale of all operators in the theory. If 
we consider the effect on Eq. (8) of a change in p 
we derive that | 


ə ə es 
(ue +B A+ 14+ 10-20) C(x ,g)=0, (9) 
where £ is given by Eq. I(4.9), and y4 (yg, Yo) is 


the anomalous dimension of the operator A (B,O), 
Thus, for example, 


(10) 


ə 
va(g*) = Me Za ina? 
Za is the renormalization constant of the operator 
A, where we have assumed that A is multiplica- 
tively renormalizable (this is true for all the op- 
erators that control the light-cone behavior of 
current products). 

In the case of interest this implies that the func- 
tions C” (q?/u?, g) satisfy 


(ug + Ble) E-ra) Č (q2/u2,2)=0, (11) 


where again we assume a unique twist-two oper- 
ator appears in the Wilson expansion of currents, 
and C™ is related, by an equation similar to Eq. 
(6), to the coefficient of the spin-n twist-two oper- 
ator. Note that the anomalous dimension of con- 
served or partially conserved currents (by which 
we mean that the dimension of the divergence of 
the current is less than four) vanishes for all g, 

The solution of this renormalization group equa- 
tion is expressed in terms of the effective coupling 
constant g(t, g) defined by 


ee 
WTO, AO, 8)=8 (12) 
where 
t=4In(-q?/p?) . (13) 
In fact 


C™ (q?/u?, g)=C™ (1, Z(t, 2)) 
t 
xexp(- f dey" (ax, 8))) » (14) 
0 


The large-q? behavior of C” will thus be deter- 
mined by the large-t behavior of at, g), which 
in turn is determined by the fixed points of the 
renormalization group Eq. (12). In an asymp- 
totically free gauge theory 


g’, g) ~ batt +b ba 3t Int +0 (1/£?) , 
t —> œ 


(15) 


where b, is given by Eq. I(4.9). The anomalous 
dimension y” will behave for small g like 


y”(a)=y"g°+0(g$). (16) 


Therefore in an asymptotically free theory 
C™ (q?/u?, g) ~ const[In(-g?)]- 
t —> œ 


x[Č® (1, 0) + O(1/in(-¢°))], 
(17) 


where 


met a 
m2,” 


a (18) 

From this expression we see that the approach 
to the asymptotic region is logarithmic, i.e., the 
corrections to the asymptotic form are suppressed 
by In(-q”). The asymptotic value of C) does not 
exhibit naive (Bjorken) scaling. Instead there are 
logarithmic deviations, whose magnitude is cal- 
culable in second-order perturbation theory. 
Furthermore the tensor and quantum number 
structure of the operator-product expansion will 
be that of free-field theory, up to logarithmic 
corrections, since C” is evaluated at g=0 on the 
right-hand side of Eq. (17). 

In general there will appear in the Wilson ex- 
pansion many operators O(™ with the same quantum 
numbers and twist. In that case a given O'”) is not 
multiplicatively renormalizable, rather one must 
take linear combinations of these to obtain oper- 
ators which are multiplicatively renormalizable 
and have definite dimensions. 

Consider the Wilson expansion 


Ja (x) Jy (- x) => Cy, (x? — 1€X, 8) 
i,k 
XOP (O)p e- oupa tirtextnt..., 
(19) 


where again we suppress the tensor indices of the 
currents. The label z runs over the set of oper- 
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ators with spin n and dimension (n +2), and the 
label k denotes the tensorial and SU(3) xSU(3) 
structure of the expansion. [In other words k 
stands for the labels a,b,1,2,3 in Eq. (3).] Since 
the operators Of can mix, a change in u will, in 
general, entail a mixing of the various oper- 
ators. The Wilson coefficients will obey the re- 
normalization group equation 


(u GEA ch} = p2 vila Ch. (20) 


The anomalous dimension is now a matrix 


(y™ (D rN 


ð 
a i) 
(z H au Lo Š 
and Z, is the matrix of wave- Rinses renormal- 
ization constants of the operators om s Lë; 


(21) 


£o »A fixed 


form of Cf} is again most easily expressed in 
terms of Z(t, g) 


ane) |r exp - [xr ase ds | \ 
J ; 


xČ® (1, Z(t, g)), (23) 


where T refers to the fact that the exponential is 
to be t ordered. In an asymptotically free theory 
Z? vanishes for large t like 5, '¢7!. The anoma- 
lous-dimension matrix y” P similarly vanishes 
like 1/t. We define 


© (g)=y™ g? +r™ (8), (24) 


where r™ (g) is of order fs for small g. Then 
the large- -t behavior of C®? "x iS given by 


H (q2/u2, 2) > a 2% | exe(= Int y®\M| 


OP. ey ORAS) ac 22 S 
f renormalized 2 foare (Zo) 122) xē® (1,0), (25) 
The solution of Eq. (20) for the Fourier trans- where 
= Int z Ini 
mM” = Texp |- f dt exp (24) 7" B(t, g))exp (- r »)| ` (26) 
o 0 0 


Thus we see that the large-q? behavior of the 
Wilson coefficients is determined by the eigen- 
values of y™. These can be calculated by eval- 
uating the wave-function renormalization matrix 
Z to second order in perturbation theory. The 
structure of Eq. (25) is more transparent if we 
write the matrix y” in terms of its eigenvalues 


y™ 


y™ => yPP, (27) 
i 


where P* are projection matrices 


PY P!= 6441 , 
De 
4 


Then we have 


(28) 


ČRa?) ~ Dime 
x PO (P'M™),; ČR (1,0), 
J 
(29) 


1 
m= 0) 
a 2b, YU. (30) 


The operator-product expansion has the form (in 
momentum space) 


=a (") 
JaJa™ > [In(—q?)]-27 >, OP [P'M™],, 
l i,j 


xČ® (1,0). (31) 


The dominant operator for large q? will thus be 
picked out by the smallest eigenvalue of y™. The 
coefficient of this leading operator is not, in gen- 
eral, determined. Although the functions C$} (1, 0) 
are known (equal to their free-field-theory values) 
the matrix M “) is not. It depends on the (unknown) 


behavior of y for large coupling constant. The only 


case in which the coefficient of the dominant operator 
is known is when the matrix P'M“” vanishes identical- 
ly for allg. This willbe the case if we are considering 
the coefficient of a conserved or partially con- 
served current or the energy-momentum tensor. 

In that case the (appropriately projected) matrix 

M equals unity and the coefficient of the current, 

or energy-momentum tensor is determined. In 
addition a”) vanishes for these operators. These 
are the only operators for which one is interested 
in the numerical value of the coefficients since 
only these operators have known hadronic matrix 
elements. Also if a particular tensorial or SU(3) 
xSU(3) structure holds for the operator -product 
expansion when g=0, it will also be valid when 
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—=Q°— œ. This is because the tensorial and SU(3) 
xSU(3) structure is totally contained in the labels 
k of Cf evaluated at g=0. 

Consequently all current-algebra sum rules 
(Adler, Gross—Llewellyn Smith, etc.) and relations 


between moments of the structure functions derived 


in the parton or naive light-cone models will be 
true asymptotic theorems in our theories. 


If]. DEEP-INELASTIC SCATTERING 
IN ASYMPTOTICALLY FREE 
GAUGE THEORIES 


We now supply the explicit results one obtains 
by applying the techniques of Sec. II to an asymp- 
totically free gauge theory. In order to ensure 
sufficient generality we will work with the La- 
grangian 


Lo=-27r Fy, FY" +BY -m)p ’ (32) 


where V,,=2,+ig0° B; is the covariant derivative, 
the o° being the matrices of the representation of 
the Lie algebra of the gauge group G in the space 
of fermions. We retain the normalization of Eqs. 
1(2.3)—(2.4) and omit for brevity the gauge-fixing 
and ghost terms needed to define (32) properly. 
For simplicity we assume that G is simple, but 
allow the fermions to transform according to an 
arbitrary representation of G. 

The class of theories described by (32) includes 
the case where G=SU(3) and y transforms as a 
direct sum of three triplets. In this case we may 
interpret G as the gauge group of “color” and the 
group which mixes the three G-triplets as ordinary 
SU(3), so the fermions are just the usual “colored” 
quarks. As was discussed in paper I, there is 
some reason to hope that this theory provides a 
model of hadrons. We will adopt it in the following 
for illustrative purposes. If m=0 we also have 
ordinary chiral SU(3). Mass terms, and in general 
symmetry-breaking terms with operator dimen- 
sion less than four, have no influence on the lead- 
ing-order effects we discuss in most of this paper. 
In particular, all symmetries broken by such 
terms will be reinstated asymptotically, so that 
for instance all the operators in a given chiral 
SU(3) xSU(3) multiplet have the same anomalous 
dimensions. 

More complicated models, with charm quarks or 
a different gauge group (so long as it is non-Abe- 
lian), are clearly allowed. It might be interesting 
to see what the effects of Higgs scalars on the 
light-cone behavior of asymptotically free theories 
are (asymptotically free theories with scalar par- 
ticles were constructed in paper I). One would be 
faced with some complications due to the presence 
of at least two dimensionless coupling constants 


and more involved mixings of lowest-twist oper- 
ators. We have not carried out this study. 

Throughout this paper we shall assume that the 
gauge group G commutes with the ordinary sym- 
metry group H of the hadronic currents, e.g., 
chiral SU(3) XSU(3). This excludes models of the 
Bars-Halpern-Yoshimura type.!5 There are good 
reasons for this assumption. The spin content of 
the fundamental charged constituents of the nu- 
cleon as measured by the Callan-Gross® sum rule 
points to spin-3 constituents. A “shielding” mech- 
anism (like that contemplated in paper I) is im- 
possible if one wishes to identify these particles 
with the observed vector or axial-vector octets, 
since the shielding mechanism yields only sin- 
glets of the gauge group for physical states. The 
alternative complicated system of Higgs scalars 
required in Ref. 15 almost certainly destroys as- 
ymptotic freedom. 

We are interested in the light-cone behavior 
of the product of weak (by this we mean electro- 
magnetic or truly weak) currents. The operators 
of twist two appearing in the operator -product 
expansion of two weak currents near the light cone 
will be 


noy .n~2 Pe a 
Ones eu, Zt STr Fp a Vu, Yena E ln 


~ trace terms, (33) 


Hye * Up zi SY Yy, Vu,’ "Va, (1 ty,) p 


-trace terms, (34) 


. 


ua pg 2E SP Ya, W, Vyp LY) 2A Y 
— trace terms (35) 
3 


where recall that V, is the covariant derivative, 
which is 0,,+7g0°B), acting on fermions and ð, 
+igT’ Bi, acting on vectors, while (34) and (35) 
represent the SU(3)xSU(3) singlet and octet pieces 
of the fermion operators. The A‘, a=1,...,8, are 
the standard à matrices of Gell-Mann. The prod- 
uct of two octet currents will in general also have 
decimet and 27-plet components; these correspond 
to higher-twist operators and are therefore sup- 
pressed, [The vertices associated with the oper- 
ators (33)—(35) are given in Fig. 1.] 

In addition to the operators 70" ,"0* (we often 
suppress in the following tensor and H indices) 
there are composite operators formed from 
Faddeev-Popov ghosts which may have the same 
twist and therefore are expected to mix with H 
singlet operators. Arguments to exclude the 
ghosts on the basis of gauge invariance are un- 
convincing because the ghosts do in fact mix with 
the above operators in off-shell matrix elements, 


| 2 


gur(Ak) +k? A, AAKT? -ku A, (Ak) Aktak 


H v 
(a) 
n 
k gabn 5) 
(b) 
n 
(c) 
n 
n-2 j j -2-j t 
-a8 Ak Ze) (Ap (A po) ¢ 428) 
b H 
Ppa Pob py. a 


all momenta flow into vertices 


FIG. 1. (a) The vertex for "OY (order 0). (b) The 
vertex for "OF (order 0). (c) The vertex for ”OY (order 
g). (d) The vertex for ”OF (order g). 


as is shown by explicit calculation. 

It is argued in detail in Appendix A that one can 
solve this problem by going to a gauge in which 
no ghosts are present, so that the mixing does not 
appear. Moreover, the anomalous dimensions of 
the gauge-invariant operators "O" ,"O" and their 
mixing are gauge-independent, so they may be 
calculated in the standard Fermi-type gauges. 

In the remainder of the text we shall take the re- 
sults of Appendix A for granted and ignore the 
ghost mixings. 

According to the prescriptions of Sec. II the 
light-cone behavior of the coefficients in the 
Wilson expansion, and thereby the high-q’ be- 
havior of the moments of the structure functions, 
will be controlled by the renormalization group 
Eqs. (20). The only missing ingredient in these 
equations is the matrix y. 

As was briefly mentioned in Sec. II, operators 
of the same structure (tensor and symmetry prop- 
erties) and twist will mix, and only certain linear 
combinations of them will be multiplicatively re- 
normalizable in the usual sense. This phenom- 
enon appears already in Fig. 2(a), where we see 
that "OY, which in lowest order vanishes between 
fermion states, acquires (logarithmically diver- 
gent) contributions in higher order. In order to 
compute the y matrix of Eq. (21) to second order 
in g we must according to Eq. (22) express the 


renormalized operators in terms of the bare ones. 
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(a) 


g 4 + crossed + Pew 


(b) 


A + P\ woe 


(d) 


FIG. 2. Graphs for computing y (radiative corrections 
to matrix elements). (a) The matrix element of "OY 
between fermion states. (b) The matrix element of 
"OF between vector states. (c) The matrix element of 
"oY between vector states. (d) The matrix element of 
"OF between fermion states. 


This is most easily done by taking matrix elements 
between vector and fermion states, since the bare 
operators have trivial matrix elements. We are 
thus lead to computing the logarithmically diver- 
gent pieces of the matrix elements of "O" and "OF 
with the same tensor structure (twist) as "Of... 
and "Of e, between these states. From this we 
determine the wave-function renormalization ma- 
trix and thereby y. The computations involve 
evaluating the graphs of Figs. 2(a)—2(d). Because 
we want to amputate the external propagators 
there will also be an order-g’ contribution to the 


. diagonal elements of y due to the possibility of 


self-energy insertions on the propagators, as in 
Eq. I(5.18). . 

To summarize the results we introduce the ma- 
trix 


tire Prk 
ny = ? (36) 


F n. V 
Yvv Yvv 


where "yf, is given by (8/8 yu) In(Z)F, and so 
forth. 

The results of the computations are (for n=even, 
n 2 2) 
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av g? TE: ENNE es 
yvr =g Oe n(n=1) (n+1) (n+2) 


+45) st (a) |, (37) 


"YFF -£ Coe -y 3 >) | , (38) 


—g? 4(n*+n+2) 


nF 

Yvy = gr? n(n+1) (n+2) T(R), (39) 
- —g? 2(n?+n+2 

vee r A?) CR), (40) 


n(n? - 1) 


where, as in paper I, C,(G) is the quadratic 
Casimir operator of G evaluated on the adjoint 
representation C,(R) is the evaluation of the quad- 
ratic Casimir operator of G on the irreducible 
representation R of G to which the fermions be- 
long,'’ and T(R) is the trace of the square of a 
matrix in the Lie algebra of the representation R. 
It should be remarked that the mixing terms (39) 
and (40) vanish for fermion operators which are 
not singlets under H (we have been suppressing 

H labels on fermion operators, they should always 
be understood). The remarkable feature of Eqs. 
(37)—(38), characteristic of gauge theories,” is the 
2/3 1/j~Iin(n) term. It arises from the graphs of 
2(c) and 2(d) with three lines coming out of the © 
vertex, which in turn corresponds to the fact that 
A, occurs with 2, in the covariant derivative. 

These values of "y satisfy the constraints of 
positivity (the smallest eigenvalue of "y must in- 
crease with n), gauge-invariance, and nonrenor- 
malization of the energy-momentum tensor, which 
provides a check on the calculation (and on the 
argument of Appendix A). 

Having calculated the y matrix we can now infer 
the properties of deep-inelastic scattering accord- 
ing to the methods outlined in Sec. II. First let us 
review the general features. 

1, The approach to the asymptotic region is 
logarithmic. In other words the leading correc- 
tions to the asymptotic forms for the moments of 
the structure functions will be suppressed by pow- 
ers of In(—q7). These corrections arise from two 
sources. First, the fact that the effective coupling 
Z*(t, g) vanishes logarithmically for large t (7?~1/t) 
means that the order-g’ corrections to the Wilson 
coefficients will be suppressed by 1/I1n(— q°) for 
large —q*. It is hard to estimate how rapidly 
these corrections vanish in a realistic model, since 
this will depend on the unknown scale (u), the 
actual value of the physical coupling constant, and 
the large-g* behavior of §(g”). The large-g? 
behavior of £ is, of course, totally unknown. If, 
for example, B(Z%’)=dg?/dt were to be linear in 


g’ for large Z°? then the effective coupling would 
decrease rapidly (like a power of — q°) from its 
physical value to small values of Z? and only then 
approach zero logarithmically. In that case one 
might understand the rapid onset of scaling. 
Additional logarithmic corrections occur when 
more than one operator can contribute to a given 
moment. This occurs in our theories for the H 
Singlet component of the structure functions. In 
that case, as explained in Sec. II, the leading as- 
ymptotic behavior of the moments will be given 
by the lowest eigenvalue of the y matrix. However, 
the other eigenvectors of y will also contribute 
terms, which will be suppressed by some (in gen- 
eral noninteger) power of In(—q’). Thus a generic 


_ structure function F(x, qg*) will satisfy 


1 
f dx x" F(x, q?) ~ 
-i 0 q 


as 


$ n- gyei. 
i 


xFi(n(- 4°) t>, 
(41) 


where the aj are proportional to the eigenvalues 

of y, and the F’s approach constants (at an unknown 
rate) as q?~-—~, ee 
_ 2, Bjorken scaling is violated by finite powers of 
logarithms. These logarithmic violations are 
readily calculated in terms of the y matrix eval- 
uated in second-order perturbation theory. Here 
we must distinguish between the strong symmetry 
group singlet and nonsinglet pieces of the structure 
functions. The latter are easy to analyze since 
there is only one operator, that given in Eq. (35), 
that contributes. The relevant anomalous dimen- 
sion is given in Eq. (38). It is, of course, common 
to all the SU(3) xSU(3) operators that appear in the 
operator-product expansion. Therefore if F(x, q?) 
stands for the nonsinglet piece of F,, xF,, or xF, 
(independent of the quantum numbers of the cur- 
rents) we have that 


2 
f axs E= @) m Rng, (42) 
0 —> 0 


where 


-SOR fasaa L 
An =220,(G) - 8T(F) (1 are) ee 


and 
Q, =-=? . 


In the “red, white, and blue” quark model [H 
=SU(3)], we have 


C,(G) = 3, 
C,(R) =$, (44) 
T(R) = 
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The numerical value of these coefficients is small, 
i.e., ANS=i8, ANS=0.36, etc. An excellent inter- 
polation formula, accurate to 1% already for n=2, 
is 


ANS œ 0.296 In(n) — 0.051 . (45) 


In treating the singlet piece of the structure func- 
tions we must take into account the mixing of the 
vector-meson and fermionic operators. This mix- 
ing occurs for H-singlet operators of both normal 
and abnormal parity. In the case of abnormal 
parity the fermion operator mixes with a vector- 
meson operator given by Eq. (33) with one Fav re- 
placed by its dual. These operators contribute 
to the parity-violating structure function F,. The 
appropriate y matrix in this case will be treated 
in a subsequent publication. Here we will only 
deal with H-singlet normal parity operators. Con- 
sequently we must diagonalize the y matrix eval- 
uated above. For n=2 this matrix has the form 

(2)_& : 


.- 2 C, (R) a : 


X “gr? T(R) T(R) (46) 


Thus the smallest eigenvalue is zero, so that the 
corresponding moments will in fact scale. If 
F*(x,q") stands for the singlet pieces of F, or xF,, 
then 


1 
J dex FQ) ~ CPE, (41 
0 —> 00 


where 


1 
An= ag A" Vee tyy 
0 


—[Cy ep ~ yey) tary es yee] 7} . 
(48) 


Since the off-diagonal matrix elements of y vanish 
rapidly (like 1/n), we have to a very good approx- 
imation (1% for n= 4) that 


Az=0, 


ao), n>2. a 


S ANS _ 
fn tn o( im 


The fact that ASAS is a simple consequence of 
positivity; however, we note that the difference be- 
tween these coefficients vanishes rapidly as n in- 
creases. Already for n=4 it is less than 1%, so 
that one can hardly differentiate between the single 
and nonsinglet parts of the structure functions 
(except of course in the behavior of the first mo- 
ment). 7 

In asymptotically free gauge theories there are 
then two stages in the approach to the asymptotic 
region. In the first stage, when we have reached 


what might be called the approximate scaling re- 
gion, the effective coupling constant becomes small 
and the one-loop approximation to the renormal- 
ization-group equation becomes valid. In this re- 
gion we have scaling up to finite powers of loga- 
rithms as previously computed. The rate of ap- 
proach to this region is not determined by the 
methods of this paper. In the second stage, which 
might be called the true asymptotic region, quan- 
tities suppressed by powers of logarithms become 
effectively zero, and the sum rules (energy-mo- 
mentum and parton-model sum rules) mentioned 
below become valid. 

3. Sum rules and relations between moments of 
the structure functions which follow from the ten- 
sorial and SU(3) xSU(3) structure of the free-quark 
model are true asymptotic theorems in our theories 
[if H=SU(3)]. This transpires because we have 
chosen the strong gauge group G to commute with 
H [say SU(3)]. Therefore the vector mesons are 
neutral with respect to the SU(3) XSU(3) charges, 
and the coefficients C® vanish when g =0. Thus 
the tensorial and SU(3) xSU(3) structure of the 
Wilson expansion for large —q? will be identical to 
that of free-field theory. 

The Adler sum rulet? is of course valid for all 
q?. The Gross—Llewellyn Smith sum rule” holds, 
since for this moment the appropriate anomalous 
dimension vanishes for all g. It, however, is 
approached logarithmically, i.e., 


f} ds F (x, g) +F! (x, a] =-6 F(a?) , 

(50) 
Fa’) pw 140 oa) 

The Llewellyn Smith relation?’ holds for individual 

moments 


f dx x" {6[F 2 (x, q°)- F?" (x, q°)] 


— x[ F3 (x, q7) - F3" (x, q7)]} 


as well as the various inequalities between mo- 
ments of the structure functions discussed by 
Llewellyn Smith?’ and Nachtmann.”! 

All these are relations or sum rules for the 
moments of the structure functions. As will be 
shown below it does not necessarily follow that the 
relations are valid for the structure functions 
themselves. Thus the quantity in parentheses in 
Eq. (51) does not necessarily vanish for a given 
x like 1/In(- q?) as q? —- œ. 

Similarly the Callan-Gross relation holds for 
moments of the structure function. In other words 
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the ratio of the moments of F= F, -—2xF, and F, 
vanishes logarithmically for large - q?: 


Jor dx x" F(x, q°) ( 1 z (52) 


ann ERER SESER ATES Ld ae 
fot dxx" F,(x,q?) @+-0 \In(-q’) 


However this does not necessarily imply that 

F(x, g?)/F,(x, 4°) vanishes like 1/1n(-q?) for fixed 
x. 
Finally we can derive, in our models, a new sum 


rule which is related to the matrix element of the 
energy-momentum tensor.” The singlet piece of 
the Wilson expansion will contain the energy-mo- 
mentum tensor 6,, with vanishing anomalous di- 
mension. According to the discussion in Sec. II 
this will appear, for large -q°, in the form 


(02, of?) pP® i r) , (53) 


where we have set P® M = P®) (since the anom- 
alous dimension of the energy-momentum tensor 
vanishes), C, is the coefficient of the fermionic 
part of 0,, (OY), and the coefficient of the vector 
contribution to 6,, (Of?) vanishes when g=0. The 
projection matrix P‘) is given from Eq. (46) as 


pon l GE ay (54) 
£ “9C,(R)+T(R) \ T(R) 2C,(R) 

so that the coefficient of 6,, = 69? + 0)?’ is given, 

for large —q’, by 


T(R) 


2C,(R) + T(R) Cote) aa 


The net effect of the mixing of operators is to 
multiply the free-field theory value by y, where 


T(R) 


”*2C,(R) + T(R) 


= for the “red, white, and blue” quark model.. 
(56) 
The contribution of the energy-momentum tensor 


to the commutator of two vector (or axial-vector) 
currents is thus (for a quark model) 
i), JI AD] Tr A) S o'f? 
X[x*Ouy (0) xy + xO, Xp 
-guv Xx Oy]. (57) 


As a consequence the singlet piece of F, for 
electroproduction will satisfy 


Lx 
f dx F3, anger (%, 9°) — (Qr, (58) 
o q 


-> me ORD 


where (Q?) is the average quark charge. In the 


“red, white, and blue” quark model this has the 
value #. This sum rule holds for any hadronic 
target. However, the corrections to it, arising 
from the nonsinglet operators, might be very 
large. They vanish for infinite q? but at a rate 
governed by ANS =0.2 (in the “red, white, and blue” 
quark model). 

Similarly in the case of neutrino or antineutrino 
the singlet contribution to F”? satisfies (setting 
the Cabibbo angle equal to zero) in the quark model 

1 
f RI REA (x, q?) Fee ir. (59) 
0 q? —> =o 
To improve the convergence one can take linear 
combinations of structure functions to obtain a 
pure SU(3) singlet: 


1 
f a(R P+Fe)-(FP +E] — $r. (60) 
0 q 


Foa 


The value of this sum rule for “red, white, and 
blue” quarks is 0.48, whereas experimentally 
one has 0.72 +0.28. 

We have determined the large-(—g*) behavior of the 
moments of the structure functions. What can one 
say about the q? behavior of the functions them- 
selves? It is useful to construct an example of a 
function with roughly the correct anomalous di- 
mensions. Consider the nonsinglet piece of 
FNS (x, q?) for electroproduction. It satisfies Eq. 
(42). If the constants C® were known then one 
could construct FMS (x, q?). Since these constants 
are not known there exist many functions which 
satisfy Eq. (42). Let us approximate A by its 
asymptotic form, ANS = æ In(z) - 8. For the “red, 
white, and blue” quark model a=0.296, B=0.051. 
Then we have to find a function FNS(x, q?) which 
satisfies 


1 
lim f dx x" F(x, q°) =c™ et BL (yn +2)7o%, 


a?—> -o “9 
L=I1n In(-q°). (61) 


In addition we shall impose Regge behavior, in the 
sense that for fixed q? we demand that F(x, q?) 
approach a constant as x— 0. 

A solution of this is provided by 


@ bt e~Inx 
Fada) dye yt, (62) 
where c is an arbitrary constant. This solves Eq. 
(61) with 
echt) 
“(n +1)(n+2) ` 


One can easily construct additional solutions by 
multiplying Eq. (62) by a polynomial in x. This 


Cc” 
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function has Regge behavior, i.e., 
0. q? ere "d -y „aL 
F(0, 4 )= aLi f ye” y 


~ [ln(- ¢?)]®. (63) 


qs 
At fixed x #0,1 it behaves, for large —q’, like 


(ols) )* ght 


al (64) 


F(x, q’) ~ ( 

Accordingly F(x, q°) will (a) for any x#0,1, 
vanish faster than any power of In(—q’) for large 
enough —q? (in fact F(x, q7)<[In(- q?) rr 2), 
and will (b) increase for intermediate values of 
- q’ for sufficiently small x and decrease for xclose 
to unity. Thetransition point, at which F does not 
change as —q? increases, is roughly x= 1/ 
[In(— q?)]?, where p is some positive constant. 
(c) Finally, F(x, q?) will show largest deviations 
from scaling in the vicinity of x=1. 

Such behavior is more general than this par- 
ticular example. Specifically (a) is a simple con- 
sequence of positivity, 


1 X+€ 
RP) dy Fly, 47) 
XmE 
sef dy y” F(y, q?) 
2elx = Ee) J, i 


S Finga), (65) 
faai 
and the fact that A, ~Inz as n= %. 

This example also illustrates the rather small 
variation in q? of the asymptotic form of the struc- 
ture functions. To test for deviations from scaling 
one will require large variations of q? and mea- 
surements in the vicinity of x=1. One might 
expect a 50% variation for x~ 0.9 as —q? increases 
from 10 to 50 BeV?. 

This example further illustrates that two dif- 
ferent solutions of Eq. (61) might have quite dif- 
ferent q? behavior for a given value of x.?? Con- 
sider the solution F’ to Eq. (61) which is given by 
Eq. (62) with c replaced by c’<c. Then, for large 


-q’, 


F(x, q’) a 2 mole), ( Sane)", 
F(x, q°) ~ [In(- g”)] ae en eee 


By making the ratio c/c’ large we can arrange for 
F/F' to increase like a large power of In(— q?) for 
x near unity. 

As a consequence the parton-model relations 
for moments of the structure functions, for ex- 
ample Eq. (51), do not imply that these relations 
are satisfied for the structure functions them- 
selves. Also the fact that the moments of the lon- 


gitudinal structure functions decrease like 

1/In(- q?) relative to the moments of the trans- 
verse structure function [Eq. (52)| does not imply 
that R(q?, x)=F,(q’, x)/F,(q?, x) decreases loga- 
rithmically for all x. Indeed it might very well 
increase for x close to unity. Without additional 
input the only reliable predictions one can make 
are with respect to the q? behavior of the moments. 


IV. CONCLUSIONS 


A crucial test of asymptotically free gauge the- 
ories of the strong interactions is the verification 
of the q? behavior of the moments of the structure 
functions. In the best of all worlds one could con- 
front an infinite number of these moments with 
the predicted asymptotic forms (which are deter- 
mined solely by the gauge group and the fermion 
representation). An additional test is provided by 
the energy-momentum sum rule derived above. 

In reality, of course, it will be very hard to 
determine the q? behavior of the moments. More 
than likely the most practical place to look for 
violations of Bjorken scaling is in the vicinity of 
threshold (x~ 1). There we expect to see the 
structure functions decreasing like ever-increasing 
powers of In(- q?). 

Whether the picture described in the paper is 
consistent with experiment is an open question. 
The fact that scaling appears to have set in already 
for rather small values of q? is neither explained 
nor contradicted by our theories. The rate of 
approach to asymptopia must be determined by 
nonperturbative methods. This problem, as well 
as the understanding of the low-energy and on- 
mass-shell behavior of the theory, requires the 
development of new theoretical techniques. 


APPENDIX A: PROBLEM OF GHOST MIXING 


One is tempted to argue that only manifestly 
gauge-invariant operators appear in the Wilson 
expansion of the product of two gauge-invariant 
operators. However, inthe case of composite oper- 
ators constructed from Faddeev -Popov ghost fields 
it is not clear how to formulate the criterion of 
gauge invariance. It does not seem that the ghost 
fields have any simple transformation properties _ 
under gauge transformations (in particular, one 
can choose gauges in which the ghosts are entirely 
absent). Moreover we will show in a simple ex- 
ample that operators constructed from ghost fields 
“mix” in a highly nontrivial way with the manifest- 
ly gauge-invariant operators of lowest twist which 
control the Bjorken limit, in the sense that the y 
matrix appearing in the renormalization- group 
equations for these operators has nonvanishing off- 
diagonal entries. This appendix is devoted first to 
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showing that a nontrivial problem is involved here, 
and second to describing that argument we used to 
convince ourselves that the calculations given in 
the main test do indeed give the correct asymp- 
totic behavior for products of currents in the 
Bjorken limit. We argue that the anomalous di- 
mensions (renormalization constants) of manifest- 
ly gauge-invariant operators are independent of the 
gauge in which they are computed, so since there 
exist gauges in which no Faddeev-Popov ghosts are 
present, the correct values are obtained in gauges 
with ghosts by ignoring the ghost mixing. 

An instructive example of a gauge field theory is 
quantum electrodynamics with the unusual gauge 
choice** 


Ə ,A —20gA, A" =0, (A1) 


where g is the coupling constant and @ is a gauge 
parameter. The Feynman rules for this theory are 
given in Ref. 24. Notice the presence of Faddeev- 
Popov ghosts. | 

For simplicity we consider the theory without 
fermions. In this case we are of course just deal- 
ing with free-field theory. The manifestly gauge- 
invariant operators of lowest twist are 

CO Poca =SF yo bu,” is eee te (A2) 

where S denotes symmetrization with respect to 
the “’s. To compute yy, the anomalous dimension 
of the vector field, to order g? we need only con- 
sider the logarithmically divergent parts of the 


two self-energy graphs with vector and ghost loops. 


The result is?’ 


Yy=0. (A3) 


We can now compute "y},, the diagonal term in 
the y matrix for the operator YO™, from the 
radiative corrections to the vertex shown in Fig. 
3(a). The result is 


Jur (Ak) +k A AAKT? -k A, (AK) Ak, (Ak) 


H o 


(a) 


(b) 


FIG. 3. Graphs for computing y}y in quantum electro- 
dynamics in the gauge of Eq. (A1). (a) The vertex for 
Voln)  (b) The radiative correction to the matrix element 
of Yo" between vector states. 


Yyy =0. (A4) 


To see whether these operators mix with the ghost 
operators we need only check that the anomalous 
dimension for "O'”, sandwiched between ghost 
states y{¢, does not vanish. A calculation exactly 
like the one just described for the vectors involv- 
ing the diagrams of Fig. 4(a) gives 


_arg2 -1 a?g? 


LaO V eaa O 
Ye -= 1672 ’ Yee ~nin—1) 167 e (A5) 


Notice for œ =0 we are in ordinary Landau gauge, 
and of course no ghost mixing can occur (there are 
no ghosts). 

This shows there is an operator that mixes with 
YO™ which to lowest order (g°) is given by (p is 
the ghost field) 


COR.. eua %) = Sp*(x)y °° 3, p(x)+O(Z). (A6) 


In addition there might be order-g terms in which 
one of the derivatives is replaced by gA,, arising 
from the diagrams in Fig. 4(b). 

To complete the picture we notice that y$,, the 
anomalous dimension of the ghost operator sand- 
wiched between vector states, does not vanish. 
This is readily seen by noticing that only Fig. 5(a) 
can give a term of the structure g „, InA? (notations 
are as in Sec. III). So we know that if we failed to 
consider other operators the y matrix 


"Yyy "YTV 
n., = 
(SE mas) a 

could not have zero eigenvalues ("yyy =0, "yy #0, 
"vic #0, sodety#0). However, we know that we 
are dealing with free-field theory, so there must 
be an operator with zero anomalous dimensions. 
Evidently there are even more operators which 
mix, and there is a complicated cancellation mech- 
anism among many operators which are not indi- 
vidually gauge-invariant. 

To verify that the same phenomenon occurs in 


(b) 


FIG. 4. Mixing of vector and ghost operators. (a) 
The zero-order contributions to 0, (b) The order-g 
contributions to °0™) , 
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Yang-Mills theory we will work in a gauge where 
there are no Faddeev-Popov ghosts.” Our gauge 
is specified by 


nat =0, (A8) 


where n is a fixed vector and a is the group index. 
(This gauge condition is not Lorentz-invariant. ) 
The absence of Faddeev-Popov ghosts is a conse- 
quence of the fact that, under a gauge transforma- 
tion 


6A, (x)= x(x) -g x), Ay@)], 


6(n"A‘,) is independent of A for fields satisfying (A8). 
In order to implement gauge conditions like (A8) 
we add a term —(1/2a)(n"A,) to the Lagrangian. 
This leads to the standard Feynman rules for Yang- 
Mills theories but without the ghosts and with the 
propagator modified to be 
py , ap? y nip’ +p'n” 


16 ab 
Dgo e- -g "” Pare Pp perpin), 
(A9) 


The theory in this gauge is not renormalizable by 
power counting if +40. From now on we will al- 
ways assume that @=0 and simplify further by tak- 
ing n” lightlike, “?=0. Then the propagator is 


nN Py +Pyny 
Omg yy eee ): 


The calculation of yy and "yọy is straightforward, 
following the pattern of Sec. III. The same dia- 
grams [Fig. 2(c)] are involved, except that there 
is no ghost-loop contribution to the self-energy. 
The results are 


Yy =—-(8?7/877)C,(G)¥, 


20 ap 


Dæ (P)= (A10) 


"Yyy = (£?/807)C,(G) 


1 4 4 | 
Cai ea n). 
(A12) 


In the calculations we have taken nk =nA =0 (k is 
the external momentum). It turns out that the 
dangerous-looking denominators in (A10) all can- 
cel from internal propagators, at least in our cal- 
culation of logarithmically divergent pieces. 
Although yy is certainly gauge-dependent, the 
nyk, are not.2” In fact the "y%, have the same val- 
ues in the ghost-free gauges as in the Fermi-type 
gauges. In the ghost-free gauges the operator- 
product expansion for gauge-invariant operators 
takes a simple form. In particular there appear 
no ghost operators. The anomalous dimensions of 


(A11) | 


ee? 


\ + ae 


A 


n 
\ 
và + crossed 
k ‘\ 
H v 
(b) 


FIG. 5. Contributions to"y$,. The coefficients of 
yy arise from diagram (a). 


‘the operators appearing in it can, however, be 


computed in any convenient gauge. This demon- 
strates why the calculations given in the text are 
relevant for determining asymptotic behavior in 
the Bjorken limit, despite the apparent problem 
of ghost mixing. The same apparent problem 
arises even in quantum electrodynamics in the 
unusual gauge (A1), but here, of course, we know 
the correct answer, and it is correctly given by 
the arguments of this appendix. 


APPENDIX B: SAMPLE CALCULATION, "yE p 


In this appendix we will give details of the calcu- 
lation of "vip. The techniques used here should 
enable the reader to duplicate without too much 
difficulty most of the numerical results in this 


paper. The calculation will also clarify the origin 
of the f}; 1/j terms in the anomalous dimensions 


which are characteristic of gauge theories and give 
pronounced violations of canonical scaling as n—%, 

The key tool in these calculations is the angular 
average integral 


2 


~(n+2yl D Emoto E min-tan? 
pairs : 


(B1) 


where on the left-hand side we have the average 
over the unit sphere in Euclidean 4-space, while 
on the right-hand side the sum runs over all pos- 
sible ways of grouping the u’s into pairs. This 
formula holds for even n; for odd n the left-hand 
side vanishes. The formula is readily proved by 
induction. 

Let us evaluate the first graph of Fig. 2(d). (A 
is an arbitrary vector which is contracted into all 
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free indices of the fermion tensor.) In the 
Feynman gauge it is, after some simple Dirac 
algebra, 


P 1 
108) | Bay FTI 


x[2p7h -4(p -A)AI(A +p)". (B2) 


We are interested in the (logarithmically divergent) 
coefficient of (A -k)""!X in the expansion of this in- 
tegral. Throwing away other terms, we write 


C,(n-1)IA, +e 


Now rotating into Euclidean space and using (B1) 
gives 


“ae Ter zC,(R (eas). (B5) 


The second graph of Fig. 2(d) gives 


ig*C,(R) fz (nyt 5 oeae > 2 D (a pa aye 
(B6) 


The reader who has worked through the previous 
example carefully should be able to do this integral 
“by inspection”; the answer is 


-2 ae 
C= iE CRA) 5 (B7) 


with the obvious definition of c,. 
The third graph gives the same number as the 
second. 


i dsp ma-i, 
Ae,- A int’ = ig?C R) f Goa pare Pe, Pe 


symbolically 
c (A +k)" Int? = ig?C,(R) 


dip 2p k-Ap AB, no 
Qn Prp- APP 
(B3) 


where A is an ultraviolet cutoff. Differentiating 
(B3) on both sides n-— 1 times with respect to k 
and keeping only logarithmically divergent parts 
gives 


_ (A PRPA - 4(p + a)s). (B4) 


Putting these results together gives 


n oe 
vir- wia- na 15) 


and with our old result Yr [Eq. I (4. 1) we get the 
announced result for "yf p. 

In this calculation it was clear that the origin of 
the terms in "y growing logarithmically with n 
arose from the replacement of ordinary by co- 
variant derivitives. In terms of the “hand-waving 
argument” of Ref. 2, it appears that covariant de- 
rivatives do not “separate” in space the fields 
they are sandwiched between. 

The method outlined here is a very powerful one 
for calculations of anomalous dimensions at the 
one-loop level. Other methods have also been pro- 
posed for this purpose.”® It is doubtful that any 
method of comparable simplicity exists for higher- 
order calculations. 
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The attractive (positive) and repulsive (negative) linear eigenvalues of the Lippmann- 
Schwinger equation can become arbitrarily small in absolute value for sufficiently strong 
interactions. In this paper we show that this cannot be the case for the Faddeev equations. 
Three-body eigenvalues arising from the repulsive part of the two-body interaction must 
have a certain minimum absolute value or the Faddeev equations can predict physical bound 
states held together by this part of the interaction. It is shown that this limitation upon the 
eigenvalues arises from the suppression in the T matrix of the repulsive part of the two- 


body potential. 


The linear eigenvalues of the Faddeev equations 
have been of considerable interest recently due 
to their importance in iterative methods for de- 
termining the three-body binding energy? and in 
studies of the three-body D function.’ In this 
paper we study the nature of the linear-eigen- 
value spectrum of the Faddeev equations in the 
presence of two-body interactions which contain 
a repulsive part. Our primary result is that, 
because the repulsion in the two-body potential 
is strongly suppressed in the two-body T matrix, 


the Faddeev eigenvalues arising from the repul- 
sion in the interaction have a certain minimum 
absolute value. Physical arguments also dictate 
this behavior, for otherwise the Faddeev equa- 
tions would predict bound states held together by 
the repulsive part of the interaction. 


The eigenvalue spectrum 


For three identical particles the linear-eigen-~ 
value problem can be expressed as? 


